Superfiuid transition in a Bose gas with correlated disorder 
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The superfiuid transition of a three-dimensional gas of hard-sphere bosons in a disordered medium 
is studied using quantum Monte Carlo methods. Simulations are performed in continuous space 
both in the canonical and in the grand-canonical ensemble. At fixed density we calculate the shift 
of the transition temperature as a function of the disorder strength, while at fixed temperature 
we determine both the critical chemical potential and the critical density separating normal and 
superfiuid phases. In the regime of strong disorder the normal phase extends up to large values of 
the degeneracy parameter and the critical chemical potential exhibits a linear dependence in the 
intensity of the random potential. The role of interactions and disorder correlations is also discussed. 
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The interplay between superfluidity, interactions and 
disorder in quantum degenerate Bose systems (the so- 
called dirty boson problem) is a central topic in con- 
densed matter physics, many aspects of which are still 
unsolved and under scrutiny. Since the seminal work by 
Fisher et al. the general understanding is that in- 
teractions are essential to stabilize the system and that 
superfluidity is lost for strong enough disorder, leading 
to a normal phase which at low temperatures is identi- 
fied with the Bose glass phase. However, a quantitative 
description in terms of the relevant parameters of the 
random potential and other matters, such as the criti- 
cal behavior and the role of dimensionality, are still open 
issues. 

On the experimental side a large body of work was 
devoted to 4 He adsorbed in porous media, such as Vy- 
cor glass and aerogels @, 0| ■ These studies investigated 
the behavior of the heat capacity and of the superfiuid 
response as well as the dynamic structure factor Q 
as a function of temperature and filling. However, no 
clear evidence was observed of a compressible Bose glass 
phase. More recently the dirty boson problem has been 
addressed using ultracold atoms, which offer unprece- 
dented control and tunability of the disorder parameters 
and of the interaction strength. Interaction effects were 
studied in disordered optical potentials even though 
the main effort has been given so far to the suppression 
of diffusion for non-interacting particles (Anderson local- 
ization) ||. 

Many relevant theoretical contributions are based on 
quantum Monte Carlo simulations of the Bose-Hubbard 
Hamiltonian with disorder [fj. In this lattice model the 
physical scenario is more involved than in continuous 
space because of the role played by commensurability 
and of the existence of the interaction driven phase tran- 
sition to the Mott insulating state. Other theoretical 
approaches make use of mean-field approximations 0, @| 




FIG. 1: (color online). Typical shape of the speckle potential 
Vdi s , with averaged value Vo = h 2 /mi 2 ., shown in the direction 
(0,0,1) of the simulation box. We also show schematically 
the value of the critical chemical potential /i c . Inset: radial 
dependence (in units of the inverse momentum cut-off A) of 
the disorder spatial autocorrelation function F. The solid 
(black) line refers to an average over many realizations of 
the random field, the (green) symbols correspond to a single 
realization. 



and are not reliable in the regime of strong disorder. 

In this Letter we report on a path-integral Monte Carlo 
(PIMC) study of an interacting Bose gas in the presence 
of correlated disorder produced by 3D optical speckles. 
This random potential is relevant for experiments and al- 
lows for an independent tuning of intensity and correla- 
tion length. By increasing the disorder strength, we find 
a sizable reduction of the superfiuid transition tempera- 
ture and the shift is larger for weaker interactions. We 
map out the normal to superfiuid phase diagram, both in 
the chemical potential vs. disorder and in the density vs. 
disorder plane. For strong disorder and in the presence 
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of small but finite interactions, the critical chemical po- 
tential varies linearly with the disorder intensity and is 
essentially independent of temperature and interaction 
strength, in agreement with the existence of a mobil- 
ity edge separating localized from extended states. In 
this regime and for chemical potentials below the critical 
value, the equilibrium state is a highly degenerate nor- 
mal gas which is expected to correspond to the Bose glass 
phase. We consider a system of N identical particles of 
mass m subject to the random field Vdis and interacting 
with a short-range pairwise potential. The Hamiltonian 
is given by: 
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The interatomic potential is modeled by a hard-sphere 
interaction: V(r) = +00 if r < a and zero otherwise, 
where the hard-sphere diameter a corresponds to the s- 
wave scattering length. The system is in a cubic box 
of volume V = L 3 with periodic boundary conditions. 
Disorder is modeled by an isotropic 3D speckle potential 
defined as follows [9(: 



V dis (r) = V 



i J dk^k)VF(k)e* k 



(2) 



where Vq is a positive constant and y>(k) = 
J dr(p(r)e~' lk ' r is the Fourier transform of the complex 
field tp(v), whose real and imaginary part are indepen- 
dent random variables sampled from a gaussian distri- 
bution with zero mean and unit variance. The function 
W(k) is a low-wavevector filter defined as: PV(k) = 1 
if k < 7rA and zero otherwise. The random potential in 
Eq. ||3J| is positive definite and the probability distribu- 
tion of its intensities is given by the normalized exponen- 



tial law P(Vdis) = e 



-Vdi./V 



/Vq. If the volume V is large 



enough the disorder Vdis is expected to be self-averaging, 
i.e. spatial averages coincide with averages over different 
realizations, and one has Vb = (Vdis) = 1/V J drVdis(r). 
The mean square displacement is also determined by the 
same energy scale: Vq = \J (V^ 2 is ) — (Vdis) 2 - The correla- 
tion length £ c is defined from the spatial autocorrelation 
function, T(r') = (Vdi s (r)Vd is (r + r')) - (Vdis) 2 , as the 
length scale for which T(£ c /2) = T(0)/2. We find the 
following relation between the correlation length and the 
wave- vector cutoff A: £ c = 1.1 /A. The length scale £ c is 
typically ~ 100 times larger than the hard-sphere diam- 
eter a, allowing for a wide range of disorder intensities 
where interaction effects are well described by the s-wave 
scattering length and the details of the interatomic po- 
tential are irrelevant. The typical box size used in the 
simulations ranges from L ~ 20£ c to L ~ 50£ c . An in- 
dication of self-averaging of disorder for these values of 
L is provided by the inset of Fig. [IJ where we show the 
comparison between the autocorrelation function T aver- 
aged over many realizations of the random potential and 



the one corresponding to a single realization. The typical 
shape of the speckle potential Vdis is also shown in Fig.Q] 
typical wells have size £ c and depth Vq. We notice that 
standard experimental realizations of optical speckles are 
2D, i.e. the speckle pattern lies in the plane perpendic- 
ular to the propagation of the laser beam. We consider 
instead a 3D pattern, having the same correlation length 
in the three spatial directions. 

The energy h 2 /m£ 2 , associated with the correlation 
length £ c , and Vb provide the two relevant energy scales 
for the disorder potential. In particular, if Vb 3> H 2 /m£ 2 
the random potential is classical in nature, with typical 
wells that are deep enough to sustain many single-particle 
bound states. The opposite regime, Vb <C H 2 /m£ 2 , corre- 
sponds instead to quantum disorder, where typical wells 
of size £ c do not have bound states and these can be sup- 
ported only by rare wells of size much larger than £ c or 
with depth much larger than Vq. 

The outcomes of PIMC simulations consist of unbi- 
ased estimates of thermal averages of physical quantities, 
using the many-particle configurations R = (ri,...,rjv) 
sampled from a probability distribution proportional to 
the density matrix p(R,R,T) = (R|e~ H / fcBT |R) at the 
temperature T . In the present study we are interested 
in the superfluid density p s , obtained from the winding 
number estimator [Toj . and in the one-body density ma- 
trix (OBDM) 711(7"), whose long-range behavior defines 
the condensate density no — lim r ^oo ni(r). Our simula- 
tions are based on the worm algorithm [111 ] , which allows 
for an efficient sampling of permutation cycles, and on the 
pair-product decomposition which is well suited for stud- 
ies of dilute systems [H}. We perform calculations both 
in the canonical (at fixed density n) and in the grand- 
canonical ensemble (at fixed chemical potential (J,) [ill ]. 

We are now in a position to discuss our results. First 
we discuss the simulations carried out at fixed density. 
The scattering length and the disorder correlation length 
are also kept fixed and for the latter we choose the value 
n£ 3 = 0.24, such that there is typically one particle in 
each small sphere of radius £ c : nAir£ 3 /2> ~ 1. Results 
for the transition temperature as a function of disor- 
der strength are shown in Fig. [2] for two values of the 
gas parameter na 3 . The transition temperature T c is 
expressed in units of T c ° = (27r/i 2 /m/c B )[n/C(3/2)] 2 / 3 , 
the critical temperature of the non-interacting gas with 
£(3/2) ~ 2.612, and the results in the absence of disor- 
der are taken from Ref. [l2T| . At na 3 = 10 -4 , there is no 
appreciable change for Vq < 1 compared to T c in clean 
systems. For larger intensities we find a sizable shift that 
is well described by a linear dependence in Vq. For a 
given strength Vb the reduction of the transition temper- 
ature is enhanced for smaller values of the gas parameter, 
consistently with the instability of the ideal Bose gas in 
the presence of disorder [3] ■ The value of T c is extracted 
from the results of the superfluid fraction p s / ' p (p = mn 
is the total mass density) , corresponding to systems with 
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FIG. 2: (color online). Superfluid transition temperature as 
a function of the disorder strength for two values of the gas 
parameter no? . Open and solid symbols refer respectively to 
T c determined from the superfluid and from the condensate 
fraction. The dashed line is the prediction of Ref. [j| at na 3 = 
10~ 4 shifted by (T c — T°) in the absence of disorder. Inset: 
scaling behavior of the superfluid density for different system 
sizes and different realizations of disorder. 



FIG. 4: (color online). Critical density as a function of the 
disorder strength for different values of temperature (in units 
of h 2 /ml 2 ) and scattering length. The horizontal arrows in- 
dicate the critical value n° of the non-interacting gas. Inset: 
Density dependence of p 3 /p [(pink) squares] and no/n [(blue) 
circles] for the values Vo = and Vb = 6 Ah 2 /ml 2 of the dis- 
order strength. Here k B T = 0.13/i 2 /mil and a/4 = 0.016. 
The vertical arrow indicates the corresponding value of the 
degenerate density n° c . 





-1 

-2 

-3 

-4 

-5 

-6 
O 

i -7 

-8 
-9 



E 



O 
> 



* k B T = 0.65 all c =1 .6 x1 0' 2 
k B T = 0.13 a/^ c =7.5x10' 3 

• k B T = 0.1 3 al£ c =1 .6 x1 0' 2 




6 8 

V [b 2 lme 2 c ] 



10 



12 



FIG. 3: (color online). Critical chemical potential (shifted by 
Vb) as a function of the disorder strength for different values of 
temperature (in units of h 2 /ml 2 ) and scattering length. The 
grey shaded area denotes the superfluid phase. Inset: spatial 
dependence of the OBDM for two values of the chemical po- 
tential slightly below and above p c . Here ksT — 0.13ft 2 /ml 2 
and a/lc = 0.016. Two different system sizes are used to 
check the role of finite-size effects. 



different particle number N, using the scaling ansatz 

JV 1/3 p.(t, N)/p = /(tN 1 ^) = /(0) + f'itytN 1 '^ + ... . 

(3) 

Here, t = (T — T c )/T c is the reduced temperature, v is 
the critical exponent of the correlation length £(t) ~ t~ v ' , 
and f(x) is a universal analytic function, which allows 
for a linear expansion around x = 0. The validity of the 
scaling behavior © is presented in the inset of Fig. [21 



where the effect of different realizations of the random 
potential is also shown. The quantity N^ 1+7, '' 3 no/n, 
involving the condensate fraction hq/u and the corre- 
lation function critical exponent r\ = 0.038 of the XY- 
model universality class, is also expected to obey a scal- 
ing relation of the form ([3]). For all reported disorder 
strengths Vo, the extracted value of the critical exponent 
v is compatible with the result v = 0.67 corresponding 
to clean systems 12|]. It is worth noting that the val- 



ues of T c , obtained from the scaling law of the superfluid 
p s /p and of the condensate fraction no/n, coincide within 
our statistical uncertainty (see Fig. [2]). In Ref. [8( the 
shift ST C = T c — T® of the superfluid transition tempera- 
ture is calculated using a perturbative approach for the 
^-correlated disorder (AVdis(r)AVdis(r')) = n6(r — r'), 
where AV dis (r) = V dis (r) - (Vd is ). The T c shift is found 
to be quadratic in k, implying for our speckle potential 
that 6T C /T? = (m 2 Vo 2 ^/^ft 4 ) 2 /[2(121og2) 3 ], where 
we used a gaussian fit to the radial dependence of the au- 
tocorrelation function T and considered the limit l c — > 0. 
We report this prediction in Fig. [2] (we also add the inter- 
action contribution not accounted for by Ref. [8( , so that 
in the clean case an exact result is reproduced) . Our data 
in the regime of very weak disorder do not have enough 
precision to allow for a quantitative comparison and di- 
verge from the theory before 5T C /T® becomes apprecia- 
ble. The effect of disorder on the critical temperature 
of a hard-sphere gas was also investigated using PIMC 
methods in Ref. [14f | where, however, no significant re- 
duction of T c was reported. For stronger intensities of 
disorder, the calculation of T c becomes increasingly dif- 
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ficult, since the dependence on the realization gets more 
important and larger systems are needed in order to have 
a satisfactory self-averaging of the random potential. 

In Fig. [3] we report results for the critical chemical po- 
tential fi c obtained from calculations carried out in the 
grand-canonical ensemble. A small change of u around 
/i c translates into a drastic change in the long-range be- 
havior of the OBDM (see inset of Fig. [3]): for n < p c 
the OBDM decays to zero and corresponds to a normal 
phase, for p > p c the OBDM reaches a constant value 
characteristic of the superfluid state. If interactions are 
small but finite, we also find that the value of p c is es- 
sentially insensitive to a change of temperature and of 
interaction strength. For weak disorder, this result is ac- 
companied by a very small critical density (see Fig. 
and corresponds to a renormalization of p c due to dis- 
order in an extremely dilute gas. For strong disorder, it 
is instead consistent with the picture of a mobility edge, 
which depends only on the parameters of the random po- 
tential and separates localized single-particle states from 
extended ones. In this latter regime we find a linear de- 
pendence of p c as a function of Vq , in agreement with the 
qualitative T = prediction of Refs. HI, 1|| in the case 
of classical disorder. 

Finally we analyze the dependence of the critical den- 
sity n c on the intensity of the random potential. The 
calculations are carried out in the canonical ensemble at 
fixed temperature and scattering length. The method 
used to determine n c is shown in the inset of Fig. |U For 
a given value of Vq one increases the density and cal- 
culates the superfluid p s /p and the condensate fraction 
no/n. The results are then fitted by a power-law de- 
pendence Ps/p ~ (n — n c ) v and no/n ~ (n — n c ) ly< - 1+ '^ 
for n > n c , where the proportionality coefficients are ex- 
pected to be non-universal parameters. In the inset of 
Fig. |4] we show the results corresponding to a configura- 
tion without disorder (Vq — 0) and with strong disorder 
(Vq = 6AK 2 /m£c). The reported values are averaged over 
a few realizations of the random potential and their scat- 
ter gives an idea of the relevance of this effect. For the 
small value of the scattering length used here, the critical 
density at Vq — coincides with the non-interacting re- 
sult n° c = C,{2>/2)(mk B T/2irh 2 y/ 2 , while for the large V 
one finds that n c is about a factor of eight greater than 
n° c . It is also worth noticing that for strong disorder one 
enters a regime where no/n is significantly larger than 
p s /p. More comprehensive results are shown in Fig. [4] 
where n c is estimated from the superfluid fraction, which 
is less sensitive to finite-size effects. The results clearly 
show an increase of the critical density as a function of 
Vo , from the non-interacting degenerate density n° up to 
values ~ 20 times larger. It is also worth noticing that 
for strong disorder, an increase of the scattering length a 
is accompanied by a decrease of n c resulting in a constant 
value of the critical chemical potential (see Fig.[3|). 

In conclusion, we have investigated the superfluid crit- 



ical behavior of an interacting Bose gas in a correlated 
random medium. In the regime of strong disorder and 
low temperatures we identify a phase, where the gas is 
both normal and highly degenerate, which should be re- 
lated to the Bose glass phase predicted at T = 0. An im- 
portant question that will be addressed in future studies 
concerns the equation of state and the thermodynamic 
properties of this exotic normal phase. 
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